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Abstract 



\ The consequences of the charge symmetry breaking effects of the mass dif- 

ference between the up and down quarks and electromagnetic effects for 
^ \ searches for strangeness form factors in parity violating electron scattering 



from the proton are investigated. The formalism necessary to identify and 



^ ' compute the relevant observables is developed by separating the Hamilto- 



nian into charge symmetry conserving and breaking terms. Using a set of 
SU(6) non-relativistic quark models, the effects of charge symmetry breaking 
Hamiltonian are considered for experimentally relevant values of the momen- 
tum transfer and found to be less than about 1%. The charge symmetry 
breaking corrections to the Bjorken sum rule are also studied and shown to 
vanish in first-order perturbation theory. 
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I. INTRODUCTION 



If one neglects the mass difference between the up and down quarks and ignores elec- 
tromagnetic effects, the QCD Lagrangian that governs hadronic physics would be invariant 
under the interchange of up and down quarks. This invariance is called charge symmetry, 
which is more restrictive than isospin symmetry which involves invariance under any rotation 
in isospin space. Small, but interesting, violations of charge symmetry have been discovered 
and are described in the reviews |T]-§1. All charge symmetry breaking effects arise from the 
mass difference between the up and down quarks and from electromagnetic effects. 

The second European Muon Collaboration EMC effect 0], the discovery that valence 
quarks carry only a small fraction of the nucleon spin, and the resulting search for strangeness 
in the nucleon has brought some attention to understanding the role of nucleonic charge 
symmetry breaking. If this symmetry holds, measurements of a parity violating electron 
left-right asymmetry in electron-proton scattering can determine new form factors whose 
origin lies only in the strange and anti-strange quarks of the nucleon |^,|^. However, the 
symmetry does not hold precisely and it is of interest to estimate how small the effects can 
be. This is especially true now that the first measurement of the proton's neutral weak 
magnetic form factor finds a value of the strange magnetic form factor that is consistent 
with zero [0]. 

Another issue concerns the momentum transfer dependence of any charge symmetry 
breaking effects. In principle, the charge symmetry breaking terms, which act as a per- 
turbing Hamiltonian, can cause the nucleon to mix with states which would otherwise be 
orthogonal. Such components could cause the form factor to have a dependence which 
could emphasize the effects of charge symmetry breaking. The purpose of this paper is to 
present arguments that such a possibility can not occur. 

It is worthwhile to discuss briefly how the assumption of charge symmetry simplifies the 
analysis of parity violating electron scattering ||^. The difference in cross section for right 
and left handed incident electrons arises from the interference of the photon and Z-boson 
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exchange terms. In particular, the photon-electron coupling is vector and the Z-electron 
coupling is axial, while the boson-proton coupling is vector. The matrix element for Z- 
boson proton coupling, Mj-i^Q"^) is given by 

M^,{Q') = (p, / I u^u - dYd I p, -^{pj I sYs I p, - 4 sin^ OwJ^jAQ')- (1-1) 

Our notation is that the | p, i) denotes a proton in an initial state with momentum and 
spin denoted by i. The terms u'y'^u and d'-f^d are evaluated at the space-time origin. The 
electromagnetic matrix element of the proton is denoted as Jp fi{Q'^), and the nucleonic term 
N = p,n is defined as 

J'N,dQ') = (^^,/ I ^uru-^drd-hrs I N,t). (1.2) 

The second term of Eq. (|1.1| ) is directly related to the strangeness of the nucleon, and is 
the new feature of parity-violating electron scattering. The third term of Eq. (|1 . 1| ) is well 
measured, but to extract the strange properties it is necessary to to determine the first term 
from independent experiments. We define this term as X'^-{Q^) with 

XfAQ') ^{PJ\ - d^d I p, i). (1.3) 

If charge symmetry holds, the (u,d) quarks in the proton are in the same wave function as 
the (d,u) quarks in the neutron, and the strange quark wave functions of the neutron and 
proton are identical. In that case 

x%m = J'.AQ') - J'nAQ'): (1-4) 

and the right hand side is well measured. We aim to study the error involved in asserting 
that the equality holds exactly. 

Here is an outline of this paper. The next section is concerned with displaying the charge 
symmetry formalism which allows a definition of the terms that cause the charge symmetry 
breaking correction to X^^{Q'^). This correction 6Xj^{Q'^) is obtained as a specific matrix 
element involving the charge symmetry breaking Hamiltonian. This formalism is general, 
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but our application involves the non-relativistic quark model. This model is well enough 
founded as to allow reasonable estimates of the charge symmetry breaking effects and is 
simple enough so that some general conclusions, that go beyond the specific calculations, can 
be drawn. Three different non-relativistic quark models are defined in Section 3. Computing 
the perturbative corrections to the form factors involves summing over all of the unperturbed 
intermediate states. This sum can be simplified by using an approximation in which the 
unperturbed Hamiltonian can be treated as a number, an average excited state mass M*, 
so that the sum over states can be performed using closure. The mass M* can be chosen 
so that the first correction to the closure approximation vanishes, with the result that M* 
depends on and on the perturbing Hamiltonian. This closure treatment is worked out in 
Section 4. The charge symmetry breaking observables are computed in Section 5. Section 

6 discusses the charge symmetry breaking correction to the Bjorken sum rule 0. Section 

7 is reserved for a summary and a discussion of the implications of the calculations. In 
addition, a comparison with other theories of nucleonic charge sjTiimetry breaking [p!0|-p!^ 
is presented. 

II. CHARGE SYMMETRY FORMALISM 

The isospin formalism is elaborated in several reviews Here we apply it to the 

nucleon and to the calculation of the quantity Xj^(Q^). The starting point is to realize the 
approximate invariance of the Lagrangian under the interchange of u and d quarks. This 
makes it worthwhile to define the charge symmetry operator, which is an isospin rotation by 
180° about the y-axis (taking the z-axis to be associated with the charge). This is defined 
by 

PlsUPcs = d, (2.1) 

with 

P,, = expi7rT2, (2.2) 
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and 

T2 = (2.3) 

where q is the hght (u,d) quark field operator. 

The Hamiltonian consists of a charge symmetry conserving term Hq and a breaking term 
Hi such that 

H^Ho + Hi, (2.4) 

with 

[Ho, Pes] ^0, (2.5) 

and 

[H, Pes] ^ [Hi, Pes]. (2.6) 

The unperturbed states are denoted by a subscript and defined by 

Ho\p,i)o^\f^^+^\p,i)o^Ei\p,i)o, (2.7) 

where M is the average of the neutron and proton masses. We work to first order in Hi 
such that the physical proton is expressed in terms of the unperturbed states by 

I P, i) =1 P, + ^ rj ^Hi I p, i)o. (2.8) 
£/j — Hq 

The quantity Aj is a projection operator on to states orthogonal to the unperturbed ground 
state isospin -doublet: 

Aj = 7- I p,i){p,i I - I n,i){n,i \ . (2.9) 

The measured electromagnetic matrix elements are then obtained using first-order per- 
turbation theory as 

J'pAQ^) =0 (P,^ I {\^l'u - \dYd){l + — 1— 2A,ifO I p,i)„ (2.10) 
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and 

JuAQ^) =0 {nj\ {^u^u - \dYd){l + -^-^2KH,) \ n,z)o. (2.11) 

We may relate the neutron and proton matrix elements by using charge symmetry which 
holds for unperturbed states: 



ra,i)o = Pes I P,«)o, (2-12) 



and which also gives 



Plu^uPcs = d^d. (2.13) 
This along with Eq. ( p.6|) allows one to obtain the relation 

PlH^Pe, = H^ + AH, (2.14) 

where 

/\H = PlHPes-H. (2.15) 

This equation is useful in identifying the charge symmetry breaking parts of the Hamiltonian 
which are relevant here. In particular, the isospin-vector operators are selected and doubled 
in taking the difference between the neutron and proton. The evaluation of /S.H will proceed 
by using the identity 

PlnPos = -f^, (2.16) 

expressed in terms of field operators 73 = / d^x {u{xyu{x) — d{xyd{x)^. In first-quantized 
notation this is: 

PLTsi{}Pcs = -Tsii). (2.17) 
Using Eqs. ( |2.1CI| ) and ( |2.11| ) and recalling the definition ( |1.3| ) of the relevant quantity 



^fiiQ"^) which involves matrix elements of the physical proton state leads to the desired 
result 



xUQ") = JvAQ") - JnAQ') + 5xi{Q% (2.18) 

where 

SX^^m =0 (P,/ I {\d^^d~\u^^u)-^2l^H I p,^)o (2.19) 
6 6 tji — UQ 

Obtaining the equation for 5X'^-[Q'^) is the main result of the present formahsm. This term 
contributes to the observed parity violating signal in just the same way as the interesting 
strangeness matrix element. It is therefore necessary to have some understanding about 
its magnitude and its dependence. The quantity 5X'^^^{Q'^) can be related to charge 
symmetry breaking modifications of the form factors Ge,m or Fi^2- In particular, the Fi so 
obtained are the same as 1/2 i^P-") of Dmitrasinovic and Pollock [g. 

III. NON-RELATIVISTIC QUARK MODELS 

The preceeding formalism is completely general. Here we adopt the view that it is 
reasonable to use a set of non-relativistic quark models to understand the rough size of 
effects at low and to make first estimates of the dependence. With these models, 
the necessary evaluations are not difficult and one gains insight into the physics of charge 
symmetry breaking. 

In non-relativistic quark models the spin and momentum proton are not related so that 
we may specify our notation by the replacement 

b,^)^|p,T), (3.1) 

for a spin up proton. The spin index will be treated implicitly so that | p, T) p)- 
The Hamiltonian is specified by a set of terms 

H = K + V,on + V,^ + Vg, (3.2) 

including the kinetic energy operator the confining potential Vcon which respects charge 
symmetry, and the residual electromagnetic Vem and gluon exchange Vg interactions. 
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We shall use Eq. ( p.l5| ) to identify the charge symmetry breaking Hamiltonian Aif as a 
sum or contributions from the different terms of the Hamiltonian. Thus we shall obtain 



AH = AK + AKm + ^Vg. (3.3) 

in which each term is obtained via the operation indicated in Eq. (|2.15| ) i.e. AK = P^^KPcs — 
K. 

Specifically, the kinetic energy term is given by 

^ = EK + ^), (3.4) 

where rrii depends on whether the i'th quark is an up or down quark. We use the notation 

Am 

rui = m+ ^-T^ii), (3.5) 

in which Am = m^ — md- Then 

K = Ko + AK, (3.6) 

with 
and 

Air = Am ^ T,{z) + ^ E 1^^3(0- (3.8) 

m 2m 

The first term of Eq. (|3.8| ) does not modify the unperturbed wave function and is henceforth 
ignored. 

The electromagnetic interaction contains charge symmetry breaking and more general 
charge dependent terms. This operator is given by 

T/ ^^r- V 2 4a(i)-a(j),. , . 

Tij 2 m'^ 6 m^ 

where Qi = ^ + |t3(«) and fij = fi — fj. The charge asymmetric part of Vem is given according 
to Eq. (|2J[5|) as 
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^yem = E(^3(^) + - ^5(r-;,)[l + ■ (3.10) 

We take the gluon exchange operator to be 

V, = -a.E A. ■ A,&(r-;,)(^ + + (3.11) 

where for three quark baryons: Xi ■ Xj = — |. The long range l/vij term respects charge 
symmetry and is not included here. Such a term is included, in principle, as part of the 
flavor independent confining interaction. The charge symmetry breaking piece of Vg is given 

by 

= "^^^ ^(^^w + + 1^(0 ■ ^mm^)- (3.12) 

O III A ^ A ^ 

We note that the short-range terms of the electromagnetic and gluon exchange operators 
are rather similar, so that we may re- write the charge symmetry breaking Hamiltonian as 

Ai/ = A/r + AVt + Al^5, (3.13) 

where 

^^^ = -?E(^3(^)+r3(J));^, (3.14) 

and 

n 2 Aw TT 2 

= (-- + 3«.-^);^ E(^3(0 + r3(j))<^(r.,)(l + 3a(^) • a(j))- (3-15) 

These expressions are used to simplify the evaluations performed in the next section. 

To proceed further we need to specify the confining potential and its ground state wave 
function. We shall use oscillator confinement for most of the calculations of this paper. Thus 

we write 

I p)o =1 (I 05) I X5)+ I 0^) I Xa)) (3.16) 



V2 



where 
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Here p = ^{ri — and A = -^(ri + r2 — 2r3), and the dependence on the position of the 
center of mass is not made exphcit. Standard [|T3| mixed symmetric spin {(ps) and isospin 
ixs) wave functions are used. The mixed anti-symmetric ones are denoted by the subscript 
A. If oscillator confinement is used, the full charge-asymmetric kinetic energy operator can 
be incorporated exactly into the operator Hq. This is the procedure of Ref. [1^. We keep 
the first-order perturbative treatment here for two reasons. First, all effects of first order in 
nT-d — can be treated in the same way; and second we wish to go beyond the effects of 
oscillator confinement. However, this difference in procedure does not lead to differences in 
the results of first order in — m„. 

The above wave function can be used to compute the electric Ge and magnetic (Gm) 
form factors. In the non-relativistic quark model these are given by the expressions: 

Ge{Q') = {P I E ^(1 + r3{i))e''-'^ I p), (3.18) 

and 

GMm = M{p I E ^(1 + r3(0)^e^^^-^- | p), (3.19) 
'I ' 

in which = <?■'?• These expressions need to be discussed because the the equations that 
relate J% f i of Eq. ( |1.2|) to the form factors Ge and Gm depend on the nucleon mass and one 
must therefore specify whether it the proton or neutron mass or the average that enters. The 
discussion in [|10| shows that in the Breit frame the quantity GmI^n is proportional to the 



matrix element of the quark magnetic moment operator ij^{i)/mi. However the discussion in 
Halzen & Martin shows that the 1/Mn factor does not appear in the definition of Gm, but 
that there is a factor of Mjy in the definition of Ge- The difference arises because of different 
choices of the normalization of nucleon spinors. One can not tell which is more appropriate 
without doing a more complete treatment in which the nucleon-spinor representation is 
derived from the quark model. Since the spin and total momentum degrees of freedom are 
uncoupled, such derivation is beyond the scope of the non-relativistic quark model. Thus 
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we simply use the average mass M in Eq. (|3.19| ). This introduces a difference between our 



approach and that of Ref . |]T0| . In principle , the differences are of order ( M„ — Mp ) / M ^ 1 . 3 x 
10^'^ and ignorable [|1^]. We remind the reader that it is the quark mass difference, not the 



neutron proton mass difference, that sets the scale of the charge symmetry breaking effects. 
The former quantity is larger than the latter because it must compensate for the effects 
of the electromagnetic interaction and the quark-mass dependence of the gluon exchange 
interaction which would cause the proton to be more massive than the neutron. The values 
that rud — triu might take in different models are discussed next. 

A. Neutron-Proton Mass Difference and Model Parameters 



The parameters of the non-relativistic quark model shall be determined from the neutron 
proton mass difference and a consideration of pionic effects. In first-order perturbation 
theory the mass difference between the neutron and the proton can be expressed as a matrix 
element of AH: 



Mn - Mp =0 {p\AH\ p)o. 



Evaluating the individual terms of Eq. (|3.13|) yields the following results: 

{rrid - niu) 



){P I AK I p)c 



(3.20) 



(3.21) 



/ I A I \ / 2 / , 5 

0{p I AVem I P)o = —5 1 - 



and 

o{p I AVg I p)o 
Adding the individual terms leads to 



Mn- Mp= {rrid - rriu) 



1 



1 



3 V 7r/5 V 12m2/3 



^3^3/2 Y^g- 



'25 



2/3m2 



^3^3/2 y 9 



a / 2 
3 V ^ 



12m2/3 



(3.22) 



(3.23) 



(3.24) 
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The parameters to be determined are (3, as, and rfi. We shall use m = 337 MeV as 
this leads to a proton magnetic moment of 2.79 n.m. The model used does not include 
pionic effects because these are essentially charge symmetric (as discussed below), but any 
consideration of the parameters should take implicit account of the pion cloud. We follow 



fluctuations converges for bag radii greater than about 0.6 fm. The importance of pionic 
effects decreases as the bag radius Rb increases. The parameter P is essentially the mean 
square radius of the nucleon (which corresponds to about 0.6 of Rb)- We use the calculation 
of the the A-nucleon mass difference as a measure of pionic effects. The gluonic contribution 
is given by: 



The physical value of this difference is taken here to be 300 MeV, but pionic effects also 
contribute. So (Ma — Mjv)^ is a fraction 7 of 300 MeV. Larger values of 7 correspond to 
smaller pionic contributions and larger values of /3. Three typical choices of parameters are 
shown in Table I. We shall investigate the charge symmetry breaking using each of the three 
models. 



We are interested in computing the charge symmetry breaking observables represented 
by Eq. ( |2.19| ). The different values of /i and the different helicities specified by the quantum 
numbers i, f can be used to specify the contributions to the electric E and magnetic M 
terms. Separating these terms and using the non-relativistic wave function allows specifies 



the ideas of the cloudy bag model [|TB| in which a perturbative treatment of pions as quantum 




(3.25) 



IV. CLOSURE APPROXIMATION 



Eq. ( pJ9D to 



SGe,m{Q'^) =0 {P I Oe,m{<i) 



A 



2/\H I p)o. 



(4.1) 







where 
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OeW) ^Y.[l-"-f) e--% (4.2) 

and 

The operator Hq removes the center of mass kinetic energy operator from Hq: 

Ho = Ho-^^^. (4.4) 

The expression depends only on internal coordinates p and A, so that the projection 
operator A does not depend on the initial and final nucleon momentum: 

A = /- I p)o o{p \ - \ n)o o{n \ . (4.5) 



The evaluation of Eq. ( ^4.1|) depends on knowing the energies and wave functions of all of 
the eigenstates of Hq. We shall replace Hq by a number M'^j^.^{Q^, AH) which is expected to 
depend on the momentum transfer, whether the electric or magnetic term is to be evaluated, 
and on the operator AH. This quantity is determined from the condition that the first 
correction to the simplification of the energy denominator by treating Hq as a number 
vanishes. This determination is accomplished by adding and subtracting M*{Q'^,AH)) to 
Hq: 

Hq = M*{Q\AH) + (Hq-M*{Q\AH)) , (4.6) 



and rewriting the energy denominator of Eq. ( p.l9|) as 



uhrQ ^ m-m\q\ah) + {m-m^Iq\ah)y^^' ~ ^^-^^ 

The requirement that the (unperturbed) ground state expectation value of the second term 
of Eq. ( ^4.7| ) vanishes leads to the result: 



Ml M\ m -M= I?. " \7^'^, \' (4.8) 





[[Oe,m{q)^Hq],AH] 


p)o 





{P 




P) 
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The use of the double commutator allows a straightforward evaluation of the various average 
masses of the excited states. Observe that the these masses depend on the operator AH 
and will be different for the different contributions to AH . It is convenient to define 

AEe,m{Q\ ^H) =M- M* ,,(Q^ AH), (4.9) 

and also to use corresponding definitions for the individual contributions to AH. 

The contributions to the electric terms can be obtained in a straightforward manner. 
One simplification is that o{p \ OE^q) \ p)o = 0. Then 

(n^\ - ' 0E{q)2AK I p)o o{p \ 0E{q)2AVL \ p)o o{p \ 0E{q)2AVs \ p)o 
''^''^^ ' AEe{Q\AK) + AEe{Q\AVl) ^ AEe{Q\AVs) ' 

and 

(n'2\ - o(P I 0M{q)A2AK \ p)o 
^^''^^ ^ - AEm{Q\AK) 
o{p I 0Miq)A2AVs I p)o o{p I 0M{q)A2AVL \ p)o 

AEm{Q\AVs) ^ AEm{Q',AVl) ■ ^ • ^ 

The evaluation of the various terms AEe^m{Q^, AH) is a straightforward but tedious 
procedure, simphfied by the feature that only the Kq part of Hq contributed to the commu- 
tator [Oe,m, Hq] Some of the relevant integrals are given in Table II. Using the AK in 
the double commutator leads to the result 

AEe,m{Q^ AK) = = -2hu. (4.12) 

mp 

That the above result must be obtained is an immediate consequence of the oscillator con- 
finement: the operator acting on the ground state leads either to the ground state or to 
the 2hLj excited state. Here the procedure of evaluating the double commutator was followed 
check on the algebraic procedure. 
The use of the long range l/vij part of the electromagnetic operator leads to the following 
result for the related average excitation energy: 

-1 |Q^e-^-^/^^^,(g^/?/2) + ^-^S.iQ^P/2) 
AEeMQ , A^l) - — e-W4^^(g2^/2)_e-W6 ' ^^.13) 
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where 



and 



1 Q 

S2{x)=4x^. (4.15) 
dx 

Note that the average excitation energy turns out to be the same for magnetic and electric 
probes. This is a consequence of the simple wave functions employed and is related to the 
feature that the electric and magnetic form factors have the same dependence. 
The low momentum transfer limit, 

hm AEeMQ',^Vl) = (4.16) 

shows that the 1/r.jj operator excites states of higher energy than does the kinetic energy 
operator. 

The using delta function contribution to AH leads to the following result: 
and the low limit is given by 

40 

The latter expression shows that the delta function operator is the most effective (of the 
ones we consider) at exciting the highest energy states. 

The terms AEem depend only on the variable Q^/?/2. If one multiplies AEe^m by fhj3 
the result is a function that is independent of the three models used here. This is shown in 
Fig. 1. Note that rfi(3 decreases by a factor of about two as one changes from model 1 to 
model 3. Thus model 1 corresponds to the smallest energy denominators. We shall display 
results for Q'^(3/2 < 10. Thus the maximum value of is 1.6, 2.2 and 3.1 GeV^/c^ for the 
models 1-3. The planned parity violation experiments are planned for values of ranging 
from about 0.1 to 3 GeVVc^ [|| . 
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V. CHARGE SYMMETRY BREAKING OBSERVABLES 



We are now in a position to evaluate the effects of charge symmetry breaking for any 
value of Q^. The charge symmetry breaking interactions AVem and AVg include two-body 
interactions that can be expected to lead to effects that fall off slowly with increasing values 
of Q^. We must compare such effects with the form factors Ge and Gm computed in the 
limit in which charge symmetry holds. This is because the gluon exchange interaction Vg 
includes a charge symmetric term which will also lead to slowly falling form factors. We 
will see that this feature of the strong form factors precludes a significant enhancement of 
charge symmetry breaking effects for even the highest values of that we consider. Thus 
the first task is to evaluate Ge,m using the wave function | p)o- 



We shall evaluate | p)o as arising from the harmonic confining potential including also 
the first-order effects of Vg. Starting with perturbation theory is reasonable because the 
first-order contribution of Vg to the nucleon mass is only -60 MeV for model 1 and -25 MeV 
for model 3. We shall see that for the range of between and 3 GeV^/c^ relevant here 
the influence of Vg on the computed form factors can be reasonably large. This is especially 
true for model 1 for which — 2.3 as shown in Table I. We find 



A. Ge,m{Q'^) With Charge Symmetry 



G'E(g') = exp (-Q2/3/6) + AG'B(g2), 



(5.1) 



with 



TT J2(g^)- J2(0)-2(J4(g^)- J4(0)) 

3m2 A^e((3^AK) 



The integrals Ji{Q'^) are tabulated in table II. 



Similarly the magnetic form factor is obtained as: 



GM(g')=/ipe-^'''/' + AGM(g'), 



(5.2) 
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where /Xp = 2.79 and 

agm{q ) = aw:at^' ^'-'^ 

with 

JsAiQ') = J3a{Q') - J3,4(0)e-«'^/^ (5.4) 

The ratios AGe,m/Goe,m, where the form factors in the absence of gluon exchange are 
given by G'os(Q^) = exp {—Q'^P/6) and Goa/(Q^) = /ipexp are shown in Figs. 2 

and 3. Both ratios vanish at = 0. Charge conservation mandates that this be so 
for the electric form factor. However, the change in the magnetic term vanishes also for 

= because of the specific simplicities in the model unperturbed wave function- the 
spatially symmetric wave function multiplies the symmetric spin-isospin wave function. The 
correction AGe is reasonably small, less than 20% for all of the values of that we consider, 
but the magnetic correction, AGm can be very large. If the absolute magnitude ratio 
is larger than about 0.3, we can expect that the perturbative treatment errs by more than 
about 10%. Hence, the largest values of ^ for which the models can be considered well 
defined are ~ 5 and 7 for models 1 and 2. We will display the charge symmetry breaking 
form factors for values ^ larger than those limits to provide information about the models, 
but the reader is cautioned against taking those results seriously. 

B. Charge Symmetry Breaking 

We are now ready to evaluate the influence of charge symmetry on the measured electric 
and magnetic form factors. We work to first order in perturbation theory (considering the 
charge symmetry conserving one gluon exchange interaction as a first order effect). The 
necessary equations ( fl.lOj ) and ( |4.11| ) are evaluated using the charge symmetry breaking 



interactions of Eqs. (^.8]), ( |3.10| ) and ( |3.12| ). The average excitation energies are given in 



Eqs. (CT),(ra and PT7| ). 
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The evaluations are straightforward, so we simply express the results. We consider the 
influence of each charge symmetry breaking interaction AK, AVem, and AVg separately. 
Thus the contribution of AK to the electric form factor is given by 

while the magnetic form factor has a term 

We see that the effects are order 5m/ AE <^ 5m/fh times a small coefficient. Furthermore, 
the dependence 6Ge,m{Q^, K) ~ Q"^ (3e~^'^^/^ is different than that of the leading order 
dominant term ~ (5e~^^^/^ and this enhances the importance of charge symmetry breaking 
at the higher values of that we consider. 

Including the effects of the electromagnetic interaction between quarks leads to the fol- 
lowing contributions to the form factors: 

^^E [^)- AE,^Q2^^Vr.) w??, AE^iQ^AV,) ^^"'^ 



and 



(5.8) 



27 [AEM^ ^Vl) ?>fn? AEm{Q\ AVs) _ 
Here negligible effects arc anticipated because of the small value of a r:^ 1/137 and because 
of the large energy denominators. These terms include the integrals J3 and J4 which fall 
much more slowly than the leading order term, recall Table II. 

Including the effects of the gluon exchange interaction between quarks leads to the fol- 
lowing contributions to the form factors: 

AG!(9)fn2^ Am 20 f^. q-^p/e _ -q^p/2a. 

[Q ) - AEeiQ^AVs)f3m^27\l7r(3^^ ^ ^' ^^'^^ 

and 



The explicit formulae show the appearance of the e"*^^^/^^ term which, at higher values of 
is much bigger than the e~^^^/^ variation of the leading order term. One might expect 
that this feature would allow the charge symmetry breaking effects to stand out. However, 
the leading order charge symmetric form factors also have a term, caused by gluon exchange, 
which also varies as e~*^^^/^^. 

The computed charge symmetry breaking electric form factors are shown in Figs. 4 and 
5 which display 5Ge/Ge as a function of (5^/2/? using G^; of Eq. Fig. 4 shows the 

three contributions to 5Ge/Ge arising, in model 1, from the individual charge symmetry 
breaking terms: kinetic energy (K) electromagnetic interaction (em) and gluon exchange 
(g). The electromagnetic term gives a negligible contribution, but the other terms can give 
contributions that are as large as 1%. The sum of the three contributions are shown in Fig. 5 
for each of the three models. The effects are largest for model 1 because of its large value 
of Us- It is possible that charge symmetry breaking could be as large as 2%. If one wishes 



to assert that only small values of as are allowed [|T9|, then the maximum charge symmetry 
breaking would be about 1%. 

The computed charge symmetry breaking magnetic form factors are shown in Figs. 6 and 
7. The ratio 5Gm/Gm is displayed as a function of Q^/2/9 where Gm is given in Eq. ( ^.21) . 
Fig. 6 shows the three contributions to SGi\j/Gm arising, in model 1, from the individual 
charge symmetry breaking terms. Once again, the electromagnetic term gives a negligible 
contribution, but the terms g and K can give contributions that are as large as 1%. In this 
case the gluon exchange and kinetic energy terms tend to cancel, with the sign difference 
arising from the different combinations of spin matrix elements appearing in the magnetic 
terms. The net result shown in Fig. 7, for each of the three models, is that the largest effects 
are less than about 1% for values of /2(3 for which the models are valid. 

It is worthwhile to examine the low effects by determining the change in the mean 
square radii caused by the different terms. The unperturbed form factors each vary as 
1 — /3(5^/6, and (3 is the mean square radius. The charge symmetry breaking terms lead 
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to behavior of the form 1 — {(3 + S(3)Q'^/6. We denote the various Sf3 according to whether 
related to the electric or magnetic terms and according to the origin of the effects. The 
results are listed in Table 111. The electric terms are much bigger than the magnetic terms, 
for which the different terms tend to cancel. Thus only 6i3e is changed in a non-negligible 
manner. For model 1, the sum of the individual contributions gives for model 1 a result 
= 0.008, which corresponds to a 1.6% change in the root mean square radius. 



C. Dependence on wave function 

The previous numerical results have been obtained using the harmonic oscillator wave 
function. Are the presently obtained very small values of the charge symmetry breaking 
effects a simple consequence of this? Another way to ask this question is: Is it possible to 
find a wave function for which the effects of charge symmetry breaking are enhanced? 

The purpose of this section is to address these questions through the use of wave functions 
other than the harmonic oscillator. Such an investigation is necessarily limited but will allow 
us to make arguments that are more general. 

We start by considering the simple wave function introduced by Henley & Miller (HM) 
First the SU(6) nature of the 3-quark wave function spin-isospin wavefunction of 



Eq. (|3.16| ) is unchanged. Then "^{p, A) is replaced by a function ^^(p^ + A^). The generaliza- 
tion is to expand the square of the wavefunction in terms of harmonic oscillator wavefunctions 
of the form given by Eq. ( |3.17|) , 



KniP, A) = / df3 gHM{f3) e-(^'+^')/^. (5.11) 



Henley & Miller chose the function g{f3) so that the resulting electric form factor of Eq. (p.l8|) 
is of the usual dipole type GEiQ"^) = jq^j;^- In this case 

gHMm{7iPf = ^A'Pexp (-A2/J/6), (5.12) 

and 
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with i? = y/ p'^ + and -ft'i(x) is a Bessel function of an imaginary argument. This wave- 
function was originally used along with a semi-relativistic Hamiltonian in which the kinetic 



plus rest mass energy is given by y/p'^ + and is not suited for calculations with the 
non-relativistic operator of eq. ( ^.4| ). This is because of the non-relativistic kinetic energy 
operator has an infinite expectation value in the wave function "^hm- This very same wave 
function was also used in Ref. ( ||2l|] ). 

We shall proceed here by using a different function g{P), one which leads to a finite 
expectation values of the kinetic energy, but which also leads to a power law falloff of the 
form factor. Using this wavefunction will allow us to see if the very small effects of charge 
symmetry are associated with the rapid Gaussian fall off of form factors obtained from the 
oscillator model. In particular we take 



which gives 



9{P) = -^e~^'^/^ (5.14) 



^'ip,\)= I d/3giP)e~^''^^'^/^, (5.15) 
Jo 



and therefore 



''W-^flAVl/jAfl). (5.16) 



The integral form (|5.15|) is useful for evaluating matrix elements of local operators such 
as expiq-fi or v{rij). For such operators the actions of taking the matrix element in a 
harmonic oscillator wave function and integrating over (3 commute, i.e. one may integrate 
the harmonic oscillator matrix element times g{l3){TTj3)^ over /3 to obtain the final answer. In 
particular, the evaluation of Eq. (|3.18|) now is the integral of exp {—Q'^P/G)g{P){7iPY which 



leads to the result 



,2 \4 
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To see how this works in studying charge symmetry breaking effects, we evaluate the 
effects of the magnetic hyperfine interaction. RecaUing Eq. ( |5.1CI| ) we now have 

<«^' - -SiAl3|^/!t? i: ^^-e---(e-'^/" - e-«"-,. (5.18, 
which is evaluated as 

1 1 



(A2 + g2/4)5/2 (A2 + g2)5/2 



(5.19) 



81 AE^{Q\ Vs) 65/2 ^3 

For large values of this form factor falls roughly as which is slower than the 
behavior of Eq. (|5.17| ). Thus it might seem that at large enough the charge symmetry 
breaking effects would dominate. This, of course, is not true. The strong form factor of 
Eq.( |5.2|) has a term AGm{Q'^) of the same momentum dependence as that of the charge 
symmetry breaking term of Eq. ([4.11| ). Thus the strong form factor would also have a 



behavior and would not be encumbered by the small factor 4?^. 

There is a general lesson that can be drawn from this exercise. Small charge symmetry 
breaking effects derived from a perturbative term in the strong Hamiltonian can not lead to 
form factors of a different asymptotic form than that of the strong form factors. 

The only possibility to get new effects is from the charge symmetry breaking in the 
kinetic energy operator AK. One might think that the kinetic energy acting on | might 
generate a state vector with different behavior. To assess the importance of AK the relevant 
expressions of Sect. IV must be re-evaluated using the wavefunction of Eq. (|5.16|) . The 
calculations are tedious but straightforward, so the results will be presented after the model 
parameters are discussed. The size parameter A is chosen so that G^; of Eq. ( |5.17|) is 
consistent with a root mean square radius of 0.83 fm. This gives A = 5.90 fm~^. One may 
compare the relevant size of our present wave function with that of the harmonic oscillator 
in another manner by computing the contribution of the kinetic energy operator to the 
neutron-proton mass difference 



2 



M,. - M, = -^^^(i^rj X'iX- ^^j . (6.20) 
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which may be equated with the harmonic oscillator result of Eq. ( p.21|) . to obtain an 
equivalent harmonic oscillator parameter f3eq in the latter such that \J~^q = 0.77 fm. The 
present wave function corresponds to a larger size than the oscillators used here. For purposes 
of estimation, we take — = 5.2 MeV, which is the value of model 1. A calculation 
of all of the relevant charge symmetry breaking terms would probably lead to a value a bit 
larger than that because a value of alphas larger than the 2.3 of model 1 would be needed 
to reproduce the A-nucleon mass splitting. 

The results are shown in Figs. 8 and 9. The computed ratios of charge symmetry breaking 
effects to charge symmetry conserving ones are displayed as 5Ge/Ge or 5Gm/Gm as a 
function of Q^/5/2 where = 0.77 fm. Here the electric and magnetic form factors have 
the functional form of Eq. ( |5.17| ). Observe that the computed ratios are once again very 
small. 



VI. BJORKEN SUM RULE 

The structure function gi{x,Q'^) can be measured in lepton-nucleon deep inelastic scat- 
tering DIS by using a polarized beam and a polarized target [Q. See the reviews Here 
X is the Bjorken variable. The dependence of gi{x,Q'^) arises from perturbative QCD 
evolution effects and from higher twist and target mass corrections. For the present purpose 
of evaluating the influence of charge symmetry breaking using non-relativistic quark models 
it is sufficient to consider the Bjorken sum rule within the framework of the naive parton 
model. 

The naive parton model interpretation of the spin-dependent DIS data is that the valence 
quarks contribute very little to the proton's spin. This startling finding motivated the studies 
of parity violation in electron-proton scattering studies discussed here. The parton model 
structure functions measure the probability for finding a quark with momentum fraction x 
in the proton and which is polarized either in the same t or the opposite [ direction to the 
proton's polarization |, and the structure functions are described by the four independent 
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parton distributions {q ± qY^x); (g ± q)^{x). The function gi{x) is given by: 

gi{x) = WelAq{x), (6.1) 
^ q 

where Ag(x) = (g^ + g^)(a^) — (^"'^ + '(l^){x) is the polarized quark distribution and Cg denotes 
the quark charge. 

In the naive parton model the integral 



Aq = dx Aq(x) 
Jo 



determines the fraction of the proton's spin which is carried by quarks (and anti-quarks) of 
flavor q. Thus, as reviewed recently p2|], one obtains 



ri,p = £ dx g,ix) = ^{pA\ (4Am + Ad + As) | p, T). (6.2) 

We shall not explicitly write the proton spin | in the following development. Thus | p) is to 
be understood as denoting | p, T). The operators Aq are the axial current operators for the 
different quark flavors, q: q'^^'^^q. The axial charge qa measured in beta decays is given by 

gA = {p\Au- Ad\p). (6.3) 

The Bjorken sum rule involves the neutron matrix element: 

Tin = —{n\AAu + Ad + As\n). (6.4) 
18 

The use of the formalism of Sect. II, and the fact that the strange quark field operator is 
not influenced by rotations in isospin space, allows one to express this quantity as a proton 
matrix element: 

Tin = 4(P I 4A(i + Am + As I p) + 4 o(p I (4Au + Ad) _ ^ ^ 2AH \ p)^. (6.5) 
Taking the difference between the equations (|6.2| ) and (|6.5|) leads to the result: 

Tip - ri„ = ^ + Ar, (6.6) 
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where 



Ar = -1 o(p I (4A7i + Ad)-^2AH \ p)o. (6.7) 
io m — xig 



The first term of Eq. ( |6.6| ) represents the Bjorken sum rule and the second term is the naive 
parton model correction to it caused by charge symmetry breaking. 

We shall use the non-relativistic quark models [^] to estimate the size of AF. In this 
case the relevant operators are: 

Au = Y.^-±^as{^, (6.8) 

i 

and 

Arf = E^^^^3(^) (6.9) 

i 

SO that 

4An + Arf = ^ E as{^) + ^ E (6-10) 

i i 

The first term does not excite the nucleon and is irrelevant. The action of the second operator 
on the nucleon leads to either a nucleon or to a A. This simplifies the calculation enormously 
since only the A intermediate state needs to be included in the sum over intermediate states 
required to evaluate Eq. ( |6.7|) . 

Then the expression for AF as obtained by using the relevant operator ( p.lOp in the 
matrix element of Eq. ( |6.71 ) is simply 



Ar = ^ o{P I E ^3(^)^3(0 '^^^^J 2 AH I p)o. (6.11) 

Only the spin dependent pieces of the charge symmetry breaking Hamiltonian AH can 
contribute to this matrix element. Thus 

-(^ - ^^)-^ E((-3(0 + nUmn,)a. . a, (6.12) 

However, the spin dot products may each be replaced by unity because 
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■ I A) =1 A) 



(6.13) 



for all pairs i,j of quarks. Furthermore the expectation value of the isospin operators 
vanishes. Consider i, j = 1, 2 and note that 



The spatial symmetry of the A and nucleon wave functions insures that this vanishing occurs 
for all pairs i, j. 

The result of this calculation is that the influence of the charge symmetry breaking 
Hamiltonian on the Bjorken sum rule vanishes. This exact is due to the use of first-order 
perturbation theory within an SU(6) symmetric wave function. 



Let's summarize. The charge symmetry breaking observables relevant for parity- violating 
electron scattering and a general formalism for their evaluation are obtained in Sect. II. 
This formahsm is just a simple way to keep track of the effects of the charge symmetry 
conserving and violating ifi Hamiltonians. The observables are evaluated using a set of 
three non-relativistic quark models, each with harmonic oscillator confinement and obeying 
SU(6) symmetry, that is defined in Sect. III. The models are distinguished by their different 
size parameters, and are required to reproduce the A— nucleon mass difference, or a size- 
dependent fraction thereof. The charge symmetric breaking effects included are the effects of 
the mass difference between the up and down quarks in the kinetic energy operator and one- 
gluon exchange interaction, and the electromagnetic interaction. One obtains a reasonable 
range of values of — needed to reproduce the observed value of the neutron-proton 
mass difference. 

The charge symmetry breaking effects are small and therefore well-treated using first- 
order perturbation theory. One must sum over an infinite set of intermediate states to carry 
out the necessary calculations. This summation is aided by the approximation of treating 




(6.14) 



VII. DISCUSSION 
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Hq appearing in the energy denominator as a constant. The relevant constant AE is chosen 
so that the first correction to the approximation vanishes; see Eq. ( |4.8| ). This means that 
AE depends on the momentum transfer and the operator that excites the proton. This use 
of a constant allows one to use closure to perform the sum over intermediate states. This 
procedure is the subject of Sect. IV. 

The evaluations are presented in Sect. V. First the electric and magnetic form factors 
of the eigenstates of Hq are obtained. The strong one-gluon exchange operator gives a high 
momentum tail which dominates the Gaussian term obtained from the oscillator wave func- 
tion. Then the influence of the three charge symmetry breaking terms in the Hamiltonian 
are evaluated for the three different models. The effects due to the u-d quark mass differ- 
ence are larger than that of the electromagnetic interaction, but are themselves very small. 
The largest of the effects we find are of the order of 1% for the change in G^; caused by 
charge symmetry breaking effects. Some larger values are shown in the figures, but these 
are for values of the momentum transfer which are outside of the regime of applicability of 
the models we use. These small values arise because of the small sizes of the basic effects: 
the ratio of the quark mass difference to constituent quark mass is about 1/70 and, the tail 
caused by the strong one gluon exchange potential makes it impossible to find a region of 
momentum transfer for which these effects can stand out. This result is not a consequence 
of the use of oscillator confinement. A different wave function, in which the square of the 
wave functions is an integral of harmonic oscillator wave functions, is also used, and very 
small effects of charge symmetry breaking are obtained. 

The charge symmetry breaking correction to the Bjorken sum rule is examined in 
Sect. VI. Here the use of SU(6) symmetric wave functions is shown to lead to a vanish- 
ing correction in first-order perturbation theory. Charge symmetry breaking therefore has 
no impact on current studies of the validity of the Bjorken sum rule. 

Next consider other computations of the effects of charge symmetry breaking on the 
nucleon. The present work is most similar to that of Ref. [|r^, and our results are consistent 
with those, except for one detail (see Sect. Ill) that depends on issues beyond the scope of 
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the non-relativistic quark model. That previous calculation is extended here by including the 
effects of the strong and electromagnetic hyperfine interaction, by studying the momentum 
transfer dependence and by using a non-oscillator proton wave function. One difference is 
that in Ref. |T0| the sum over intermediate excited states is saturated by the A(1550). We 
use closure to carry out the sum. The charge symmetry breaking operators are isovector 
which act on T=l/2 states so that the intermediate states can have either T=l/2 or T=3/2. 
The closure approximation used here allows the T=l/2 states to be included. 

Let's also discuss the work of Ma who who presents his results in the form —SG%j ~ 
0.006 0.088 n.m.. This is small compared to the current experimental error of about 
0.2 n.m., but relatively important compared to the rather small strange magnetic moment 
G^_f=-0.066 n.m. from the baryon-meson fluctuation model of Ref. p4|. The abstract states 



that the neutron proton mass difference of leads to an excess of n = 7r~p over p = n'^n 
fluctuations, but two different effects actually lead to the results. 

The first is claimed to arise from the light cone treatment of the non-interacting propaga- 
tor. If the light cone treatment is used, a term — Ai^ with Ai^ = Z^Li ^^C'"' replaces 
our non-relativistic inverse propagator M — Hq. The value of P^P = M^, the square of 
the nucleon mass. In Ma's treatment this takes on the two values of or m^. This leads to 
a numerical result r^^^ = P{p = 'K^n)/P{n = ir^p) = 0.986, which corresponds to an excess 
of 0.2% of n = 7r~p fluctuations if P{p = 7T~^n) ^ P{n = n^p) ^ 0.15. However, one should 
perform a light-cone perturbation theory treatment of the charge symmetry breaking, which 
involves treating the Hamiltonian operator P+P^ as a sum of charge symmetry breaking 
and charge symmetry conserving terms. In this case the relevant eigenvalue, analogous to 
M used here in Sect. II, must be |(mp + m^). Using this value and changes the above result 
of 0.986 to 0.992; the effect is reduced by a factor of 2. 

Actually, the biggest effect used by Ma is caused by the assumption that the radius P ~ 1 
fm of the nn'^ component of the proton is 2.5% smaller than that of the pvr" component 



of the neutron. This according to |]TT| could be caused by Coulomb effects. However the 



effects of the Coulomb potential and electromagnetic interactions in loop graphs is of order 
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a/ii <^ .025, so that Ma's effect, while physically reasonable, is estimated to have too large 
a value. A reasonable estimate of the effect is could be 0.03 n.m [^. This effect is too small 
to be relevant to experiments. 

Finally, consider the work of Celenza and Shakin [|1^] who computed the deep inelas- 
tic structure functions of the nucleon using a quark model which preserves translational 
invariance. Effects of charge symmetry breaking enter in their calculation of the ratio of 
F2{x)/F2{x). They can reproduce the experimental values of this ratio by allowing the 
neutron confinement radius to be about 10 percent larger than the corresponding proton 
radius. Such an effect is well motivated, but the 10% value is much larger than the < 1% 
effects found here. 

The net result is that the effects of charge symmetry breaking on the nucleon wave 
function can be expected to be very small. Only a limited number of models are discussed in 
the present paper, but it seems very difficult to construct a reasonable model of the nucleon 
which incorporates large charge symmetry breaking effects. That any charge symmetry 
violating effect in the Hamiltonian Hi has its analog in the symmetry preserving Hamiltonian 
Hq is a model independent statement. Thus the present conclusion about the lack of import 
of charge symmetry breaking effects seems to be true in independent of the particular model 
used. 

This work was stimulated by a talk given by D.H. Beck at the national INT. I thank 
A. Bulgac for making some very useful comments. This work is partially supported by the 
USDOE. 
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FIGURES 

FIG. 1. Energy denominators vs Q'^[5/2. The term m of the figure is fh of the text. SoUd-energy 



denominator for the long range operator Eq. (4.13). Dashed-energy denominator for the short range 
operator Eq. ( 4.17| ). 



FIG. 2. Changes in electric form factors due to Vg. The ratio of the second AGa/ to first 
//pC"*^^^/^ = Gqm terms of Eq. ( |5.2| ). The numbers refer to the models 1-3. 

FIG. 3. Changes in magnetic form factors due to Vg. The ratio of the second /S.Ge to first 
g-<9^/3/2 = Gqe terms of Eq. (^.l|). The numbers refer to the models 1-3. 

FIG. 4. Charge symmetry breaking electric form factor. The different contributions are shown 
for model 1 

FIG. 5. Charge symmetry breaking electric form factors for each of the three models. 

FIG. 6. Charge symmetry breaking magnetic form factor. The different contributions are shown 
for model 1 

FIG. 7. Charge symmetry breaking magnetic form factors for each of the three models. 
FIG. 8. Use of the wave function of Eq. (5.16). Change in electric form factor. 
FIG. 9. Use of the wave function of Eq. (5.16). Change in magnetic form factor. 
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TABLES 



TABLE I. Parameters of the Non-Relativistic Quark Models 



Model 


1 


2 


3 


V^(fm) 


0.7 


0.6 


0.5 




1.20 


1.02 


0.85 




2.3 


1.20 


0.35 


md - rriu (MeV) 


5.2 


3.8 


2.3 


7 


0.80 


0.67 


0.33 



TABLE II. Relevant Integrals 



O 



1 



ri2 

03 = ^ 

04 = 5(ri3) 



Ji(Q') 

J4(g') 



2 -Q2/3/24 

4 V 
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TABLE III. Charge symmetry breaking changes in mean square radii, Am = m„ — 



Cause 




AT/ 






—1 Am 




—as Am 1 f2 




3 m 


8 "J mv^VT 


model 1 


0.0051 


-1.14x10"^ 


0.0029 


model Z 


O.OUdo 


i.o7 xiU 


0.0014 


mouei o 


u.uuzz 




U.UUUO 


Cause 




AT/ 


AT/ 


Mm. 


9 m ^ 


27 V T 


-7as Am 1 /2" 
240 m T 


model 1 


0.0024 


-0.0013 


-0.00069 


model 2 


0.0013 


-0.0011 


-0.00031 


model 3 


0.0005 


-0.0004 


-0.00007 
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